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Dissipative particle dynamics (DPD) is a novel particle method for mesoscale modeling of complex 
fluids. DPD particles are often thought to represent packets of real atoms, and the physical scale 
probed in DPD models are determined by the mapping of DPD variables to the corresponding 
physical quantities. However, the non-uniqueness of such mapping has led to difficulties in setting 
up simulations to mimic real systems and in interpreting results. For modeling transport phenomena 
where thermal fluctuations are important (e.g., fluctuating hydrodynamics), an area particularly 
suited for DPD method, we propose that DPD fluid particles should be viewed as only 1) to provide 
a medium in which the momentum and energy are transferred according to the hydrodynamic laws 
and 2) to provide objects immersed in the DPD fluids the proper random "kicks" such that these 
objects exhibit correct fluctuation behaviors at the macroscopic scale. We show that, in such a 
case, the choice of system temperature and mapping of DPD scales to physical scales are uniquely 
determined by the level of coarse-graining and properties of DPD fluids. We also verified that DPD 
simulation can reproduce the macroscopic effects of thermal fluctuation in particulate suspension 
by showing that the Brownian diffusion of solid particles can be computed in DPD simulations with 
good accuracy. 



Dissipative particle dynamics (DPD) is a method 
developed primarily for the simulation of complex 
fluids at mesoscopic scales 1 -^. While it is often 
thought that DPD beads represent packets of real 
atoms moving collectively, the statistical mechani- 
cal foundation of such a view for DPD model re- 
mains obscure 4 . Associated with the ambiguity of 
the exact nature of DPD beads is the ambiguity in 
mapping of DPD scales to the physical scales. While 
the mapping of length is straightforward, the map- 
ping of time is more complicated. Depending on 
the problems being studied, mapping of DPD time 
to physical time has been based on bead diffusion 
rate, bead thermal velocity, or externally imposed 
time scalea^jl, to name just a few. In all but a few 
casesi, temperature of the DPD system was cho- 
sen arbitrarily. Established procedures for mapping 
between DPD and physical scales and for choosing 
system temperature are not yet available. 

Our interest in DPD originates from the need to 
study transport phenomena in particulate suspen- 
sion where thermal fluctuations may play a crit- 
ical role&S. Such transport can be described by 
the fluctuating hydrodynamics theories 1 ^. However, 
solving the hydrodynamics equations in particulate 
suspension is computationally demanding and in- 
troducing thermal fluctuations that rigorously sat- 
isfy the Fluctuation-Dissipation Theorem (FDT) is 
challenging 11 . DPD is well-suited for studying such 
phenomena as momentum/ energy conservation and 
FDT are guaranteed by the way DPD models are de- 
signed, and colloidal particles can be modeled eas- 
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ily by bonding DPD beads together— However, 
there are two unresolved issues in DPD simulation 
of fluctuating hydrodynamics, namely, how to set up 
a model for a given physical system and how to in- 
terpret the simulation results. For example, to study 
the diffusion of a 30 nm diameter particle in water at 
300 K, what should be the temperature of the DPD 
system and how to map the results to dimensional 
values are not clear. 

To address the above issues, we propose to aban- 
don the idea that DPD fluid beads are "clumps" 
of real fluid atoms, but view them together as a 
"media" that provides an arena for the transport 
of momentum, energy and particulates that satisfies 
the fluctuating hydrodynamics laws. This idea is 
inspired by a recent paper on coarse-graining in col- 
loidal suspensions 4 . To bridge the DPD and physical 
systems, we require that 

1 . The time scale of diffusional transport of mo- 
mentum (and energy if in non-isothermal sim- 
ulations) inside the DPD fluids should match 
that of the real fluids; 

2. DPD fluids should provide objects immersed 
in them the proper random "kicks" such that 
these objects exhibit correct fluctuation be- 
haviors at the macroscopic scale. 

We now consider the application of the above re- 
quirements in the simulation of colloidal particles 
(diameter: d) dispersed in fluids (density: p, kine- 
matic viscosity: V, temperature: T). Properties of 
corresponding DPD fluids are denoted by the same 
symbol as in real fluids but without the bar, e.g., 
density of DPD fluids is denoted as p. We will limit 
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our discussion to isothermal simulations, and the ex- 
tension to non-isothermal simulations is straightfor- 
ward. The DPD model reads^ 

dr l = Vidt; mdvi = F?dt + Ffdt + Ff\fdt (1) 



where to, r^, and v 



are the mass, position, and ve- 
locity of bead i, respectively. T is the system tem- 
perature. Fp, Ff and Ff are the conservative, dis- 
sipative and random forces acting on bead i, respec- 
tively. These forces are given by^ 1 ^ 



^a ii w(r ij /r c )e i 



(2) 



F? = J2~7v w d( r ii/ r i)( e ir v ij) e ij ( 3 ) 
F? = XJ^i^(^7^)%ey- (4) 

where ay is the conservative force coefficient and 
r ij = \ r ij\ ~ \ v i ~ r j|- w anc ^ w d are weighting func- 
tions with cutoff distances of r c and r~, respectively, 
ey = /r^, and «y — Vj — Vj. #y is a symmetric 
random variable with zero mean and unit variance. 
jij and (Jij are lated by the Fluctuation-Dissipation 
Theorem as Yy — a 2 j/2k B T, where k B is the Boltz- 
mann constant. Mass, length, and time in the above 
model are measured by to, r c , and \Jk B T '/to, respec- 
tively. We assume that the unit length and time in 
DPD map to physical length [L] and time [t], re- 
spectively. To satisfy the proposed requirements, we 
enforce 



d 2 /v[t] = a f/l 



^/k B T/pd^[L\/[t] = ^Jk B T/pd 3 



Combining Equs. ((5][6|) 

[t] = 
k B T = 



pv 2 k B T 



pv 



(5) 
(6) 

(7) 
(8) 



Equs. ([7][5]) can be used to setup and analyze DPD 
simulations of fluctuating hydrodynamics. Equ. I[7| 
provides the mapping of time scale in DPD model. 
Equ. (H|) indicates that the temperature in DPD 
simulation of fluctuating hydrodynamics is deter- 
mined by the level of coarse-graining (represented 
by [L] and p) and transport properties of DPD flu- 
ids (represented by v). If v and p are known, then 
the system temperature is uniquely determined. In 
practice, since v of DPD fluids is a function of their 
temperature, the temperature of DPD system can 
only be determined after the dependence of v on 
temperature is known. While Equs. {5{ and ([6]) have 
been used to map DPD models to physical scales, 



enforcing them simultaneously in fluctuating hydro- 
dynamics simulations and thus leading to a unique 
choice of system temperature has not been proposed. 

To demonstrate the above mapping scheme and 
to investigate the ability of DPD in modeling fluc- 
tuating hydrodynamics, we study the diffusion of a 
single nanoparticle immersed in fluids with p — 10 3 
kg/m 3 and D = 0.89 xI(T 6 m 2 /s at 300 K. We set 
[L]=10 nm and p = 6.0. The nanoparticle is built 
by bonding 117 DPD beads together and is modeled 
as a rigid body. In Equ. ([2]), ttf(ry/r c ) = 1 — r.ij/r c 
with r c = 1.0 and ay is set to a// = 10.0 and af p = 
17.0, where // and fp denote fluid-fluid and fluid- 
nanoparticle interactions, respectively. In Equ. ([3]), 
w d( r ij/ r i) — \/l — Tij/r d suggested in Ref. [IH is 
used with rf = 1.0. ay in Equ. jU is taken as 5.0 
for all bead pairs. Using these parameters, we first 
measured the viscosity of DPD fluids as a function 
of temperature, and the temperature in DPD model 
of the nanoparticle-fluids system is then determined 
via Equ. §5§ to be 0.3875. To investigate the abil- 
ity of DPD model in reproducing the macroscopic 
effects of thermal fluctuations, we computed the dif- 
fusion coefficient (-Dp) of the nanoparticle, which 
is a macroscopic manifest of the thermal fluctua- 
tions. Fig. QJa) shows the mean square displace- 
ment of the particle, and a D p of (1.59±0.12)x 10~ 4 
was obtained. We also computed D p independently 
by using the Einstein-Stokes law D p = k B T /Q-k Rp,, 
where p is the fluid viscosity. From the particle- 
fluids pair correlation function shown in Fig. QJb), 
the nanoparticle radius R was determined to be 1.65, 
with which a D p of 1.55xl0 -4 was then computed. 
Given the ambiguity in the definition of particle di- 
ameter and the statistical uncertainty of DPD re- 
sults, the agreement between DPD simulation and 
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FIG. 1: (a). Mean square displacement of the nanoparti- 
cle immersed in fluids, (b) Particle-fluid pair correlation 
obtained from DPD simulation. 
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Einstein-Stokes prediction is very good. This veri- 
fies the ability of DPD in capturing the macroscopic 
effects of thermal fluctuations in particulate suspen- 
sions, which has been assumed but not explicitly 
confirmed in the literature. 

In summary, we proposed a method of mapping 
DPD simulation of fluctuating hydrodynamics, and 
the temperature in such simulations is determined 



uniquely by the level of coarse-graining and proper- 
ties of DPD fluids. Following the proposed method 
of choosing system temperature, we studied the 
Brownian diffusion of a nanoparticle and showed 
that DPD simulation can reproduce the macroscopic 
fluctuation of nanoparticle immersed in fluids with 
good accuracy. 
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